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The in-plane correlation lengths and magnetic disorder of magnetic domains in a transition metal 
multilayer have been studied using neutron scattering techniques. A new theoretical framework is 
presented connecting the observed scattering to the in-plane correlation length and the dispersion of 
the local magnetization vector about the mean macroscopic direction. The results unambiguously 
show the highly correlated nature of the antiferromagnetically coupled domain structure vertically 
throughout the multilayer. We are easily able to relate the neutron determined magnetic dispersion 
and domain correlations to magnetization and magnetotransport experiments. 
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The interplanar coupling and in-plane magnetic domains are essential to an understanding of the origin of the large 

giant magnetoresistance effect (GMR) m in magnetically coupled multilayers. This coupled with the advances in thin 

film deposition techniques O has led to a huge interest in magnetic multilayer systems specifically with respect to 

their device application possibilities. The GMR effect arises from the antiferromagnetic (AF) coupling of typically a 

transition metal ferromagnet (e.g. Co) across a noble metal non-magnetic spacer (e.g. Cu). This AF coupling can 

be realized by tuning the noble metal spacer thickness 0]. The change in resistivity results from the spin dependent 

scattering of the conduction electrons which depends not only on the magnetic moment alignment but also on the 

interfacial disorder H and the magnetic domain structure. Until recently the question of the relationship between 

fj ■ magnetic domain structure and interlayer coupling has not been explored experimentally. It is clear that a vertically 

l — ' \ incoherent magnetic domain structure will have the effect of lowering the GMR by preventing perfect AF alignment 

in adjacent layers 0. In studies of a weakly coupled system such as [Cu(60A)/Co(60A)]x20 it was shown that the 

reduction in the GMR from the as-prepared state to the coercive state can be understood as a loss of vertical coherence 

I/-) ■ of the AF coupling ||. The situation is different in the strongly coupled samples investigated here, the results of 

which clearly show magnetically correlated domains at the coercive field which extend vertically throughout the entire 

multilayer. 

The investigation of structurally rough interfaces is well established and makes use of diffuse x-ray scattering tech- 
niques. The theoretical tools for analyzing various surface morphologies are well advanced [J7|-|l0|l- Recent advances in 
x-ray techniques have applied this structural formalism to the study of magnetically rough systems |pdj-[l7|| . Neverthe- 
^^ ■ less, the problem of quantifying magnetic roughness remains difficult primarily due to the indirect and complicated 
nature of the spin-photon interaction |18|.[19| . This problem can be resolved by neutron techniques for which the direct 
C^ ■ interaction between the neutron's dipole moment and the sample magnetization is well understood. 

In this letter we have performed neutron scattering measurements on magnetically coupled multilayers and quan- 

,— ■ titatively determined the field dependence of the magnetic roughness and domain distribution. The large lateral 

coherence length of the neutron beam (> 30/im p0[) ensures that the measurements sample many magnetic domains. 

O ■ Since the neutrons are highly penetrative the measurements also sample the whole multilayer vertically unlike the 



transition metal Lm x-ray measurements 21 which sample primarily the uppermost interfaces because of the high 
x-ray absorption coefficient. 

We prepared Co/Cu and Co/Ru multilayers of 50 bilayer repeats, with Cu and Ru spacer thicknesses corresponding 
to the 1st and 2nd AF maxima of the coupling oscillation, for different thicknesses of the magnetic layer. The samples 
were deposited by dc magnetron sputtering in a custom vacuum system with a base pressure of 2xl0 -8 Torr. The 
multilayers were grown on 20 mmx25 mm pieces of (001) Si wafer with the native oxide layer left intact. The working 
gas was 3 mTorr of Ar, and deposition rates for Co, Cu and Ru were all ^3 A/s. Smaller 10 mmx2 mm samples were 
grown in the same growth run for magnetoresistance and Magneto-Optic Kerr Effect (MOKE) measurements. 

The reflectivity measurements, both polarized and non-polarized, were performed on the time-of-flight polarized 
neutron beam reflectometer CRISP at the ISIS facility, Rutherford Appleton Laboratory |22|]23J 1. To maximize the 
flux at the sample position, for the diffuse scattering measurements, the reflectometer was run in a non-polarized 
mode with an incident wavelength range of 0.5 A-6.5 A. An electromagnet at the sample position provides an in-plane 



reversible field of ±7 kOe. The scattered neutrons are detected by a 1-dimensional 3 He detector. The combination 
of the time-of- flight technique and the multidetector ensure that both the parallel (Qz) and perpendicular (Qx) (to 
the surface normal) components of the neutron wave- vector transfer (see fig flj) are obtained in a single measurement. 
Typical acquisition times are of the order of 2 hours for an entire reciprocal space map, which compares favorably 
with resonant x-ray techniques KM . 

Fig|j](a) presents the observed reciprocal space intensity map for the nominal [Co(20A)/Cu(20A)] x50 multilayer at 
remanence. This Cu thickness corresponds to the 2nd AF ordering peak. Although we have similar data for other Co, 
Cu and Ru layer thicknesses we shall concentrate on this sample in this article. Three features are apparent in the 
data: the specularly reflected ridge (Qx—0 A -1 ), the first order nuclear Bragg peak (Qz— 0.15 A -1 ) and the Bragg 
peak corresponding to the AF periodicity (Qz— 0-075 A -1 ). This peak is entirely magnetic in origin. The first order 
Bragg peak indicates that the bilayer thickness is ~42 A. The narrow width in Qz (see inset) implies that the AF 
order is coherent throughout the whole multilayer. 

A major conclusion of this paper results from the comparison of the Qx distribution of the two peaks. The nuclear 
Bragg peak is sharp but the AF peak is diffuse. The roughness is therefore predominantly magnetic. We associate 
this magnetic roughness with the existence of AF coupled domains. The diffuse scattering is strongly peaked in Qz 
therefore our data gives evidence for the coherent coupling of the magnetic domains vertically through the multilayer. 
Note, no evidence for diffuse scattering from uncorrelated regions was observed which would be uniformly distributed 
in Qz f24j. Applying a saturating field (Figj|(b)) destroys the AF correlations resulting in a ferromagnetic alignment. 
Fig.g details sections in Qx through the AF Bragg peak as one cycles from close to remanence to positive saturation 
and then reverse to negative saturation. At low fields (<100 Oe) the scattering is dominated by the diffuse scattering. 
As the field is increased to saturation only the specular ridge remains. Equivalent sections through the nuclear Bragg 
peak reveal no evidence of diffuse scattering. 

In order to quantitatively analyze our data and to characterize the domain structure of the multilayer system, we 
now present a new theoretical framework for diffuse magnetic scattering in systems with a spatially inhomogeneous 
magnetization profile, m(r), where r is the position vector. Considering a system where m(r) is in-plane, we can 
write m(r) = mo (cos </>(r), sin <j>(r), 0), with phase angle, </>(r), and amplitude, m . Thus, we consider solely directional 
variations ofm(r) which describe the different orientations of the magnetic domains. We treat m(r) and therefore 4>(y) 
as random variables, characterized by the correlation function C(|r|) = (<p(r) (f)(0)) . This is similar to the treatment 
of structurally rough surfaces by Sinha et al. H. Thus, <f>(r) plays a role reminiscent to the local height variation in 
the non-magnetic case and we parametrize C(r) as: 

C(r) = a 2 exp (-r/£) . (1) 

a = ((f) 2 ) is the width of the angular distribution and therefore characterizes the local magnetic roughness. £ is the 
lateral correlation length, i.e. a measure for a typical domain size. We consider the magnetic scattering function within 
the Born approximation, <S(Q) oc Y) a g J d 3 re lC *' r I 5 a f) — Q a Q/3 ) (flo (r)ro/3(0)), where Q a is a unit vector component 
of the transferred momentum, Q. Performing the average with respect to the different domain orientations by 
assuming a Gaussian distribution for <f>(r), we find in addition to the specular scattering, iS spec .(Q) = rn^e^ 17 <5(Qii), 
the diffusive scattering function: 

Sdiff . (Q) = m 2 e- a2 J d 2 v e lQ " r [(1 - Q 2 X ) sinh (C(r)) 

+ (l-Q|,)2sinh 2 (C(r)/2)j . (2) 

Here, Q|| is the in-plane component of Q. In our experimental geometry the detector aperture is set up such that 
the neutron intensity is integrated out over Qy, which is parallel to the applied field. Finally, when evaluated at the 
AF ordering vector, it holds that ttiq = /isin(6*/2), where fj, is the Co magnetic moment and 8 the angle between 
Co moments in adjacent layers. Before we analyze our data using Eq. we emphasize that, within the Gaussian 
approximation, we treat the angle 4>(r), normally restricted to ±7r, as an unrestricted variable. Therefore, we cannot 
describe a system with equally distributed angles, i.e. with (e 1 ^ 1 -^) = 0. However, this practically never occurs 
after the system was exposed to an external field, even if this field is set to zero or equal to the coercive field, see 
also our results in Fig. where the AF peak always has some specular component. Furthermore, our result for the 
diffusive magnetic scattering, Eq. 0, is anisotropic with respect to Q x and Q y . This is because magnetic fluctuations 
perpendicular to the field are not only larger in amplitude but also more extended in space compared to those parallel 
to the field. In the present scattering geometry this effect is negligible since Q 2 , Q 2 <C Q 2 ~ 1. However, for smaller 



Q z or by using polarized neutrons, we predict a pronounced anisotropy of the diffuse magnetic scattering if 5(Q) is 
averaged with respect to the component of Q parallel or perpendicular to the field. 

The results of numerically convoluting the specular and diffuse (Eq.0) contributions with the instrumental reso- 
lution function and performing a least-squares fit to the data are shown in Fig. 0. The agreement between theory 
and experiment is excellent. Panels (a) and (b) display the observed magnetization loop as measured by MOKE 
measurements and the normalized change in resistivity respectively. The MOKE loop and magnetoresistance curve 
both indicate good AF coupling. For fields close to rcmancncc the Co layers have a global anti-parallel alignment (c) 
with a large magnetic roughness (d) and a characteristic domain size of « 1/xm. For increasing fields three effects 
occur. The anti-parallel alignment across the non-magnetic spacer is diminished. The orientational domain distribu- 
tion within a given layer focuses around the applied field direction and the domain size increases to « 7/im. Even at 
reasonably large fields the moments are not perfectly aligned about the field direction. At +200 Oe there still remains 
a substantial domain distribution although the orientation of adjacent layers is nearly ferromagnetic (mo — ► 0). At 
these fields the diffuse scattering approaches the experimental background (primarily from incoherent scattering) and 
represents the limits of the current measurements. For this reason we cannot measure values of a close to zero in 
Fig. Q as saturation is approached. At saturation only the structural specular peak remains. The data clearly show 
the hysteresis in moving around the loop. The fact that all quantities in Fig. [| follow this hysteresis loop reveals the 
close correlation between the GMR effect and the magnetic domain correlations. The decrease of the angle 8 between 
the magnetization in neighboring layers causes not only the expected decrease of the GMR but also a decreasing 
roughness accompanied by an increase of the domain size. The latter effect is particularly drastic for fields above 
the coercive field, H c , i.e. £ w £of(H/H c ) with non-linearly growing function f(x). Below H c the dominant effect 
is the focusing of domain orientations. The observed increase in domain size from ~1.5 fim to 7 /im extracted from 
the data is typical of such systems p5pq ] . Interestingly, qualitatively similar effects have been observed around the 
nuclear/ferromagnetic Bragg peak for an equivalent ferromagnetically coupled system. 

We now turn to the domain reversal mechanism. This can take place by directional magnetization fluctuations 
which are focused by the applied field, or by domain wall motion where the axis of magnetization is unchanged due to 
a dominating magneto-crystalline anisotropy. In this latter scenario a would be a measure for the probability of the 
magnetization aligning anti-parallel to the external field rather than the width of a directional distribution. However, 
our samples are only weakly anisotropic. Imaging of 2nd AF coupled multilayers p9] strongly supports a domain 
reversal predominantly via rotation of the magnetization in agreement with our analysis. 

Our results can be compared with the recent work of Borchers et al. (Rcf. [pi) on a weakly coupled system. Their 
data show that the reduction in the GMR from the as-prepared state to the coercive state can be understood as a 
loss of coherence of the AF coupling. The more strongly coupled samples investigated here clearly show magnetically 
correlated domains at the coercive field. By increasing the exchange coupling by replacing Cu with Ru |27J there is a 
significant change in the domain structure. Even at remanence the correlation length is £ = (7 ± 3) x 10^A, with the 
diffuse tail extending across the whole of the Qx range. It is important to note that only with the very long lateral 
coherence lengths of the neutron beam used is it possible to measure £ accurately when it is so large. A full discussion 
of the dependence on spacer material and magnetic and spacer layer thicknesses will be published separately |28|1 . 

To summarize, within a new theoretical framework we have quantified the magnetic domain structure in an AF 
coupled multilayer using diffuse magnetic neutron scattering. The systematic study of the field dependence of the 
diffuse scattering reveals a close relationship between magnetic roughness, domain size, interlayer coupling and the 
GMR effect itself. 
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FIG. 1. (a): The observed scattering from the [Co(20A)/Cu(20A)] x50 multilayer in zero applied field. The intensity centered 
at Qz =0.075 A -1 corresponds to the AF ordering wave- vector and arises purely from the magnetic ordering. The intensity at 
twice this wave-vector is the first order multilayer structural Bragg peak. The dark areas represent the kinematical limits of 
the measurement, (b) The corresponding measurement in a saturation field of H=700 Oe. The AF correlations are suppressed 
leaving only the specular ridge (Qx=0) and the first order Bragg peak. The inset shows the specular reflectivity for the low 
(open symbol) and high (closed symbol) field data. 
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FIG. 2. The diffuse scattering observed at the AF peak as a function of applied field. Each scan is offset for clarity. 
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FIG. 3. (a): The room temperature MOKE magnetization loop for the [Co(20A)/Cu(20A)]x50 sample, (b) The magnetore- 
sistance. Panels (c,d,e) represent the parameters described in the text. The lines are simply guides to the eye. 



